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The noise kernel is the vacuum expectation value of the (operator-valued) stress-energy bi-tensor 
which describes the fluctuations of a quantum field in curved spacetimes. It plays the role in 
stochastic semiclassical gravity based on the Einstein-Langevin equation similar to the expectation 
value of the stress-energy tensor in semiclassical gravity based on the semiclassical Einstein equation. 
According to the stochastic gravity program, this two point function (and by extension the higher 
order correlations in a hierarchy) of the stress energy tensor possesses precious statistical mechanical 
information of quantum fields in curved spacetime and, by the self-consistency required of Einstein's 
equation, provides a probe into the coherence properties of the gravity sector (as measured by 
the higher order correlation functions of gravitons) and the quantum and/or extended nature of 
spacetime. It reflects the medium energy (referring to Planck energy as high energy) or mesoscopic 
behavior of any viable theory of quantum gravity, including string theory. The stress energy bi-tensor 
could be the starting point for a new quantum field theory constructed on spacetimes with extended 
structures. In the coincidence limit we use the method of point-separation to derive a regularized 
noise-kernel for a scalar field in general curved spacetimes. It is useful for calculating quantum 
fluctuations of fields in modern theories of structure formation and for backreaction problems in the 
early universe and black holes. One collorary of our finding is that for a massless conformal field 
the trace of the noise kernel identically vanishes. We outline how the general framework and results 
derived here can be used for the calculation of noise kernels for specific cases of physical interest 
such as the Robertson- Walker and Schwarzschild spacetimes. 



*Electronic address: Nicholas.G.Phillips@gsfc.nasa.gov 
t Electronic address: hub@physics.umd.edu 



1 



I. INTRODUCTION 



The central focus of this work is the noise kernel, which is the vacuum expectation value of the stress-energy 
bi- tensor for a quantum field in curved spacetime It plays the role in stochastic semiclassical gravity [HJ|] 

similar to the expectation value of the stress-energy tensor in semiclassical gravity [Q. We believe that this two 
point function (and the hierarchy of higher order correlation functions) of the stress energy tensor possesses precious 
statistical mechanical information about quantum fields in curved spacetime and reflects the low and medium energy 
(referring to Planck energy as high energy) behavior of any viable theory of quantum gravity, including string theory. 
The context for our investigation is stochastic gravity and the methodology is point-separation. Let us first examine 
the context and then the methodology. 



A. Stochastic Gravity 



In semiclassical gravity the classical spacetime (with metric g ab ) is driven by the expectation value () of the stress 
energy tensor T ab of a quantum field with respect to some quantum state. One main task in the 70's was to obtain 
a regularized expression for this quantum source of the semiclassical Einstein equation (SCE) (e.g., ^|,|5| and earlier 
work referred therein). In stochastic semiclassical gravity of the 90 's the additional effect of fluctuations of 
the stress energy tensor js| [To| enters which induces metric fluctuations in the classical spacetime described by the 
Einstein-Langevin equation (ELE) |]ll| , |i2| . This stochastic term measures the fluctuations of quantum sources (e.g., 
arising from the difference of particles created in neighboring histories flll l) and is intrinsically linked to the dissipation 
in the dynamics of spacetime by a fluctuation -dissipation relation p3|,p.4[ , which embodies the full backreaction effects 
of quantum fields on classical spacetime. 

The stochastic semiclassical Einstein equation, or Einstein-Langevin equation, takes on the form 

G ab [g] + A ffQb = 8^G(T ab c + T ab qs ) 

T«J^(T ab ) q + T° b (1.1) 

where G ab is the Einstein tensor associated with g ab and A, G are the cosmological and Newton constants respectively. 
Here we use the superscripts c, s, q to denote classical, stochastic and quantum respectively. The new term T* b = 2r ab 
which is of classical stochastic nature measures the fluctuations of the energy momentum tensor of the quantum field. 
To see what T ab is, first define 

t ab {x) = f ab (x) - (f ab (x))I (1.2) 

which is a tensor operator measuring the deviations from the mean of the stress energy tensor. We are interested in 
the correlation of these operators at different spacetime points. Here we focus on the stress energy (operator-valued) 
bi-tensor t ab (x)t c <d'{y) defined at nearby points {x,y). metric.) A bi-tensor is a geometric object that has support at 
two separate spacetime points. In particular, it is a rank two tensor in the tangent space at x (with unprimed indices) 
and in the tangent space at y (with primed indicies). 
The noise kernel N ab c'd' bitensor is defined as 

4AWd'0r,2/) = \ {{i ab {x)Xjd>{y)}) (1.3) 

where {} means taking the symmetric product. In the influence functional (IF) jlj] or closed-time-path (CTP) 
effective action approach jl6| the noise kernel appears in the real part of the influence action (l^Jl^]. [] The noise 



1 The derivation of noise kernel from IF or CTP coarse-grained effective action is well-known from earlier work ]ll| , [lli||ll| .p|] . 
Generally speaking the effective action approach excels in the treatment of backreaction effects of quantum fields on background 
spacetimes since it has self-consistency built in. However, most calculations of the influence actions derived so far require some 
form of perturbative expansion, because it is difficult to obtain the influence functional in closed form, (exceptions are only for 
spacetimes with high symmetry, and for special classes of coupling between the field and spacetime). We adopt the covariant 
point separation method because we can obtain a expression for the noise kernel in general curved spacetime with explicit 
dependence on the world function between two points and its covariant derivatives. This will be useful for later exploration of 
extended feature of spacetime. 
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kernel defines a real classical Gaussian stochastic symmetric tensor field r a b which is characterized to lowest order by 
the following relations, 



(Tab(x)) T = 0, {Tab{x)T c , d ,{y)) T = N abdd , (x, y) , 



(1.4) 



where ( ) T means taking a statistical average with respect to the noise distribution r (for simplicity we don't consider 
higher order correlations). Since T a b is self-adjoint, one can see that N a b c 'd' is symmetric, real, positive and semi- 
definite. Furthermore, as a consequence of (1.3) and the conservation law W a T a b — 0, this stochastic tensor r ab is 
divergenceless in the sense that V a T a b = is a deterministic zero field. Also g ab T a b(x) — 0, signifying that there is no 
stochastic correction to the trace anomaly (if T a b is traceless). (See (7|). Here all covariant derivatives are taken with 
res pect to the background metric g ab which is a solution of the semiclassical equations. Taking the statistical average 
of (1.1) with respect to the noise distribution t, as a consequence of the noise correlation relation (1.4), 



(TZb)r = (T ab ) c , 



(1.5) 



we recover the semiclassical Einstein equation which is (1.1) without the T^ b term. It is in this sense that we view 
semiclassical gravity as a mean field theory. 



B. Relation to Semiclassical and Quantum Gravity 



Stochastic semiclassical gravity thus ingrains a relation between noise in quantum fields and metric fluctuations. 
While the semiclassical regime describes the effect of a quantum matter field only through its mean value (e.g., 
vacuum expectation value), the stochastic regime includes the effect of fluctuations and correlations. We believe 
precious new information resides in the two-point functions of the stress energy tensor which may reflect the finer 
structure of spacetime at a scale when information provided by its mean value as source (semiclassical gravity) is no 
longer adequate. To appreciate this, it is perhaps instructive to examine the distinction among these three theories: 
stochastic gravity in relation to semiclassical and quantum gravity [^,^| . The following observation will also bring out 
two other related concepts of correlation (in the quantum field) and coherence (in quantum gravity). 



1. Classical, Stochastic and Quantum 



For concreteness we consider the example of gravitational perturbations h a b in a background spacetime with metric 
g a b drive n by the expectation value of the energy momentum tensor of a scalar field as well as its fluctuations 
t a b(x) p3p4| l. Let us compare the stochastic with the semiclassical and quantum equations of motion for the metric 
perturbation (weak but deterministic) field h. (This schematic representation was made by E. Verdaguer in Q). The 
semiclassical equation is given by 

Uh = lQnG(T) (1.6) 

where () denotes taking the quantum average (e.g., the vacuum expectation value) of the operator enclosed. Its 
solution can be written in the form 



h = JC(T), hih 2 = I I CiC 2 (T)(T). (1.7) 
The quantum (Heisenberg) equation 

ah = IQttGT (1.8) 

has solutions 

h = jcf, (hh) = J J dCiiff)^ (1.9) 

where the average is taken with respect to the quantum fluctuations of both the gravitational (<?) and the matter (</>) 
fields. Now for the stochastic equation, we have 

ah = 167rG((f) + r) (1.10) 
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with solutions Q 

k = J C(f)+ J CT, h 1 h 2 = J [ CiC 2 [{f){f) + ((f)T + T(f))+TT] (1.11) 

Now take the noise average () T . Recall that the noise is defined in terms of the stochastic sources r as 

(r) T =0, {T 1 T 2 ) T = {f 1 f 2 )-{f 1 ){f 2 ) (1.12) 

we get 

{hih 2 ) T = f [ C 1 C 2 {ff) i (1.13) 



Note that the correlation of the energy momentum tensor appears just like in the quantum case, but the average here 
is over noise from quantum fluctuations of the matter field alone. 

2. Fluctuations, Correlations and Coherence 

Comparing the equations above depicting the semiclassical, stochastic and quantum regimes, ^ we see first that in 
the semiclassical case, the classical metric correlations is given by the product of the vacuum expectation value of the 
energy momentum tensor whereas in the quantum case it is given by the quantum average of the correlation of metric 
(operators) with respect to the fluctuations in both the matter and the gravitational fields. In the stochastic case the 
form is closer to the quantum case except that the quantum average is replaced by the noise average, and the average 
of the energy momentum tensor is taken with respect only to the matter field. The important improvement over 
semiclassical gravity is that it now carries information on the correlation of the energy momentum tensor of the fields 
and its induced metric fluctuations. Thus stochastic gravity contains information about the correlation of fields (and 
the related phase information) which is absent in semiclassical gravity. Here we have invoked the relation between 
fluctuations and correlations, a variant form of the fluctuation-dissipation relation. This feature moves stochastic 
gravity closer than semiclassical gravity to quantum gravity in that the correlation in quantum field and geometry 
fully present in quantum gravity is partially retained in stochastic gravity, and the background geometry has a way 
to sense the correlation of the quantum fields through the noise term in the Einstein-Langevin equation, which shows 
up as metric fluctuations. 

By now we can see that 'noise' as used in this more precise language and broader context is not something one 
can arbitrarily assign or relegate, as is often done in ordinary discussions, but it has taken on a deeper meaning 
in that it embodies the contributions of the higher correlation functions in the quantum field. It holds the key to 
probing the quantum nature of spacetime in this vein. We begin our studies here with the lowest order term, i.e., the 
2 point function of the energy momentum tensor which contains the 4th order correlation of the quantum field (or 
gravitons when they are considered as matter source) Q Progress is made now on how to characterize the higher order 
correlation functions of an interacting field systematically from the Schwinger-Dyson equations in terms of 'correlation 
noise' |lj],|2(J, after the BBGKY hierarchy. This may prove useful for a correlation dynamics /stochastic semiclassical 
approach to quantum gravity ||. 

Thus noise carries information about the correlations of the quantum field. One can further link correlation in 
quantum fields to coherence in quantum gravity. This stems from the self-consistency required in the backreaction 



2 In this schematic form we have not displayed the homogeneous solution carrying the information of the (maybe random) 
initial condition. This solution will exist in general, and may even be dominant if dissipation is weak . When both the 
uncertainty in initial conditions and the stochastic noise are taken into account, the Einstein - Langevin formalism reproduces 
the exact graviton two point function, in the linearized approximation. Of course, it fails to reproduce the expectation value of 
observables which could not be written in terms of graviton occupation numbers, and in this sense it falls short of full quantum 
gravity. This comment was made by E. Calzetta to the author of [[|. 

3 The observations in this section were first made in Q 

4 Although the Feynman- Vernon way can only accomodate Gaussian noise of the matter fields and takes a simple form for 
linear coupling to the background spacetime, the notion of noise can be made more general and precise. For an example of 
more complex noise associated with more involved backreactions arising from strong or nonlocal couplings, see Johnson and 
Hu@ 
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equations for the matter and spacetime sectors. The Einstein-Langevin equation is only a partial (low energy) 
representation of the complete theory of quantum gravity and fields. There, the coherence in the geometry is related 
to the coherence in the matter field, as the complete quantum description should be given by a coherent wave function 
of the combined matter and gravity sectors. Semiclassical gravity forsakes all the coherence in the quantum gravity 
sector. Stochastic gravity captures only partial coherence in the quantum gravity sector via the correlations in 
the quantum fields. Since the degree of coherence can be measured in terms of correlations, our strategy for the 
semiclassical stochastic gravity program is to unravel the higher correlations of the matter field, starting with the 
variance of the stress energy tensor and through its linkage with gravity, retrieve whatever quantum attributes (partial 
coherence) of gravity left over from the high energy behavior above the Planck scale. Thus in this approach, focussing 
on the noise kernel and the stress energy tensor two point function is our first step beyond mean field (semiclassical 
gravity) theory towards probing the full theory of quantum gravity. 



C. Point Separation and Noise Kernel 



So far we have explained the physical motivation for investigating the noise kernel and stress energy bi-tensor. We 
now turn to the methodology. 

In the light of the above discussions, the point separation scheme introduced in the 60's by DeWitt [gl| will be well 
suited for our purpose here. It was brought to more popular use in the 70's in the context of quantum field theory in 
curved spacetimes p2|-p4[ as a means for obtaining a finite quantum stress tensor. Since the stress tensor is built from 
the product of a pair of field operators evaluated at a single point, it is not well-defined. In this scheme, one introduces 
an artificial separation of the single point x to a pair of closely separated points x and x' . The problematic terms 
involving field products such as 4>(x) 2 becomes <p(x)4>(x'), whose expectation value is well defined. One then brings 
the two points back (taking the coincidence limit) to identify the divergences present, which will then be removed 
(regularization) or moved (by renormalizing the coupling constants), thereby obtaining a well-defined, finite stress 
tensor at a single point. In this context point separation was used as a technique (many practitioners may still view 
it as a trick, even a clumpsy one) for the purpose of identifying the ultraviolet divergences. 

By contrast, in our program, point separated expression of stress energy bi-tensor have fundamental physical 
meaning as it contains information on the fluctuations and correlation of quantum fields, and by consistency with 
the gravity sector, can provide a probe into the coherent properties of quantum spacetimes. Taking this view, we 
may also gain a new perspective on ordinary quantum field theory defined on single points: The coincidence limit 
depicts the low energy limit of the full quantum theory of matter and spacetimes. Ordinary (pointwise) quantum 
field theory, classical general relativity and semiclassical gravity are the lowest levels of approximations and should 
be viewed not as fundamental, but only as effective theories. As such, even the way how the conventional point- 
defined field theory emerges from the full theory when the two points (e.g., x and y in the noise kernel) are brought 
together is interesting. For example, one can ask if there will also be a quantum to classical transition in spacetime 
accompanying the coincident limit? Certain aspects like decoherence has been investigated before (see, e.g., p5|), 
but here the non-local structure of spacetime and their impact on quantum field theory become the central issue. 
(This may also be a relevant issue in noncommutative geometry). The point- wise limit of field theory of course 
has ultraviolet divergences and requires regularization. A new viewpoint towards regularization evolved from this 
perspective of treating conventional pointwise field theory as an effective theory in the coincident limit of the point- 
separated theory of extended spacetime. This is discussed in our recent paper on fluctuations of the vacuum energy 
density and the validity of semiclassical gravity (lQl . 

The paper is organized as follows: In Section III, we review the method of point separation. In Section III, we 
discuss the procedures for dealing with the quantum stress tensor bi-operator at two separated points and the noise 
kernel. We derive a general expression for the noise kernel in terms of the quantum field's Green function and its 
covariant derivatives up to the fourth order . (The stress tensor involves up to two covariant derivatives.) This result 
holds for x 7^ y without recourse to renormalization of the Green function, showing that N a b c ifi' (x, y) is always finite 
for 1/1/ (and off the light cone for massless theories). Usingthis result, we show there is no stochastic correction 



to the trace anomaly, in agreement with results arrived at in [[11 14 1. In Section IV we briefly review how "modified" 
point separation [^6| is used to get a regularized Green function. This amounts to subtracting a locally determined 
"Hadamard ansatz" from the Green function. With this we show how to compute the fluctuations of the stress tensor, 
culminating in a general expression of the noise kernel and its coincident form for an arbitrary curved spacetime. Paper 
II in this series will apply these formulas to ultrastatic metrics including the Einstein universe (^] , hot flat space and 
optical Schwarzschild spacetimes. Paper III treats noise kernels in the Robertson- Walker universe and Schwarzschild 
black holes, from which structure formation from quantum fluctuations [j29f and backreaction of Hawking radiation 
on the black hole spacetime pCj can be studied. Finally, we intend in Paper IV to use the symbolic routine to derive 
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or check on analytic expressions for the stress-energy bi-tensor in de Sitter and anti-de Sitter spacetimes. (Martin, 
Roura and Verdaguer [£3l| have obtained analytic expressions for Minkowski and restricted cases of de Sitter spaces.) 
The former is necessary for scrutinizing primordial fluctuations in the cosmic background radiation while the latter 
is related to black hole phase transition and AdS/CFT issues in string theory. 



II. POINT SEPARATION 



We start with a short overview of point separation in this section. This paper will present the general schema 
and the results for the coincident limit of the noise kernel. The details, along with how they are carried out using 
symbolic computation, will be described in a separate paper by one of us |32] ]. Here we first review the construction 
of the stress tensor and then derive the symmetric stress tensor two point function, the noise kernel, in terms of 
the Wightman Green function. Our result is completely general for the case of a free scalar field (the separation is 
maintained throughout) and thus can be used with or without regularizing the Green function. Of course, without use 
of a regularized Green function, the y —> x limit is divergent. We also give the noise kernel for a massless conformally 
coupled free scalar field on a four dimensional manifold. We then compute the trace at both points and find this 
double trace vanishes identitically for the massless conformal case; implying that there is no noise associated with the 
trace anomaly. This result holds separate from the issue of regularization. 

Having the point-separated expression for the noise kernel and using the regularized Green function, as is done for 
the stress tensor, we obtain an expression for which the y — > x limit is meaningful. We use the so-call "modified" 
point separation prescription |26| |2Sfj . In this procedure, the naive Green function is rendered finite by assuming the 
divergences present for y — ► x are state independent and can be removed by subtraction of a Hadamard form. We 
review this prescription and use it to obtain a definition of the noise kernel for which the y — > x limit is meaningful, 



A. n-tensors and end-point expansions 



An object like the Green function G(x, y) is an example of a bi-scalar. it transforms as scalar at both points x and 
y. We can also define a bi-tensor T ai ,.. an y ...y (x,y): upon a coordinate transformation, this transforms as a rank n 
tensor at x and a rank m tensor at y. We will extend this up to a quad-tensor T ai ... ani b' 1 ---b' n2 c'{---c'^ 3 d'/' - d^ which 
has support at four points x, y, x' , y', transforming as rank m, ri2, n.3, 714 tensors at each of the four points. This also 
sets the notation we will use: unprimed indices referring to the tangent space constructed above x, single primed 
indices to y, double primed to x' and triple primed to y' . For each point, there is the covariant derivative V a at that 
point. Covariant derivatives at different points commute and the covariant derivative at, say, point x' , does not act 
on a bi-tensor defined at, say, x and y: 

Tab'-c-d' = T a b'-d';c and T a b'- C " = 0. (2.1) 

To simplify notation, henceforth we will eliminate the semicolons after the first one for multiple covariant derivatives 
at multiple points. 

Having objects defined at different points, the coincident limit is defined as evaluation "on the diagonal" , in the 
sense of the spacetime support of the function or tensor, and the usual shorthand [G(x, y)] = G(x, x) is used. This 
extends to n-tensors as 



-^ai"'a n , bi---b' c'J---c ! ' d'J 

1 "11 "2 1 "3 1 



I" ...jm 



^ai"-a ni bi---b n2 ci---c n3 d±---d n4 : (^-^) 



i.e., this becomes a rank (n\ + n 2 + n% + rii) tensor at x. The multi-variable chain rule relates covariant derivatives 
acting at different points, when we are interested in the coincident limit: 



5 We know from work in the 70's in quantum field in curved spacetime [Q-^) that there are several regularization methods 
developed for the removal of ultraviolet divergences in the stress energy tensor. Their mutual relations are known, and 
discrepancies explained. This formal structure of regularization schemes for quantum fields in curved spacetime should remain 
intact as we apply them to the regularization of the noise kernel in general curved spacetimes. Specific considerations will 
of course enter for each method. But for the methods we have employed so far, such as zeta-function, point separation, 
smeared-field flluj applied to simple cases (Casimir, Einstein, thermal fields) there is no new inconsistency or discrepancy. 
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This result is referred to as Synge's theorem in this context. (We follow Fulling's H discussion.) 

The bi-tensor of parallel transport g a b is defined such that when it acts on a vector vy at y, it parallel transports 
the vector along the geodesies connecting x and y. This allows us to add vectors and tensors defined at different 
points. We cannot directly add a vector v a at x and vector w a ' at y. But by using g a b , we can construct the sum 

v a + g a b Wb' ■ We will also need the obvious property g a b = g a b . 

The main bi-scalar we need for this work is the world function o~(x, y). This is defined as a half of the square of the 
geodesic distance between the points x and y. It satisfies the equation 



a = ~a*o- w (2.4) 



Often in the literature, a covariant derivative is implied when the world function appears with indices: a a = o~' a , 
i.e. taking the covariant derivative at x, while a a means the covariant derivative at y. This is done since the vector 
—a a is the tangent vector to the geodesic with length equal the distance between x and y. As a a records information 
about distance and direction for the two points this makes it ideal for constructing a series expansion of a bi-scalar. 
The end point expansion of a bi-scalar S(x, y) is of the form 

S(x, y) = A<® + a p A^ + a p a q Af q + o p o q a T Af qr + a p a q a r a s A^ r3 + • • • (2.5) 

where, following our convention, the expansion tensors j4l™... Qn with unprimed indices have support at x and hence the 
name end point expansion. Only the symmetric part of these tensors contribute to the expansion. For the purposes 
of multiplying series expansions it is convenient to separate the distance dependence from the direction dependence. 
This is done by introducing the unit vector p a = a a />/2a. Then the series expansion can be written 

S(x, y) = A<® + aiA^ + a A& + Ja& + a 2 A^ + ■■■ (2.6) 

The expansion scalars are related to the expansion tensors via 
A {n) = 2 n / 2 A ( pX- Pn P Pl ■■■v Pn - 

The last object we need is the VanVleck-Morette determinant D(x,y), defined as D(x,y) = — det (— a :a b')- The 
related bi-scalar 

Al/2 = / D(x,y) \ " 

satisfies the equation 

A 1 / 2 (4 - a, p p ) - 2A 1 / 2 , p a' p = (2.8) 

with the boundary condition [A 1 / 2 ] = 1. 

Further details on these objects and discussions of the definitions and properties are contained in |^3| and [j32| . 
There it is shown how the defining equations for a and A 1 / 2 are used to determine the coincident limit expression for 
the various covariant derivatives of the world function ([<r ;o ], [<J :a 6], etc.) and how the defining differential equation 

for A 1 / 2 can be used to determine the series expansion of A 1 / 2 . We show how the expansion tensors A^.., an are 
determined in terms of the coincident limits of covariant derivatives of the bi-scalar S(x,y). Ref. |||] details how 
point separation can be implemented on the computer to provide easy access to a wider range of applications involving 
higher derivatives of the curvature tensors. We will say a few words about this to end this section. 



B. Symbolic Computation 

Since the noise kernel involves up to four covariant derivatives, we expand the Green function out to fourth order 
in the distance between the points in an end-point expansion. One of the main advantages of this approach is that 
it readily lends itself to implementation in a symbolic computing environment |S4| , [33| . So the task becomes one of 
developing the necessary series expansions of the various geometric objects, including the conformal transformation 
properties. From the symbolic computational viewpoint, the key to this procedure is to recognize that all the objects 
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for which a series expansion are needed are defined by first or second order covariant differential equations. This allows 
us to develop recursive algorithms for the symbolic manipulations. The lowest order terms in the series expansions are 
known [ p3"|]35| ] . For example, Christensen's work on the stress tensor goes to second order in a a (ending up creating 
results spanning a page and a half in his manuscript). The noise kernel results get much longer. (Some expansion 
coefficients end up over 400 terms in length.) 6 

We have developed a systematic set of computer code, using MathTensor [p3[ , running in the Mathematica pifl 
environment on a workstation, to carry out much of the work contained in here and the following papers in this series. 
We use it to compute the fluctuations of the quantum stress tensor for scalar fields from a given form of the metric 
and an analytic expression for the Green function. We can also deal with cases where the Green function is only 
known in an analytic form in a conformally related spacetime. No recourse is made to numerical methods; the results 
are exact up to the accuracy of the given analytic form of the Green function before regularization. The code itself 
carries out the regularization via the "modified" point separation prescription. 



III. STRESS ENERGY BI-TENSOR OPERATOR AND NOISE KERNEL 



Even though we believe that the point-separated results are more basic in the sense that it reflects a deeper structure 
of the quantum theory of spacetime, we will nevertheless start with quantities defined at one point because they are 
more familiar in conventional quantum field theory. We will use point separation to introduce the biquantities. The 
key issue here is thus the distinction between point-defined (pt) and point-separated (bi) quantities. 

For a free classical scalar field with the action 

S[4>] =~f {m 2 <f> 2 + 4? RH + 4>-, v #*) V9d A x. (3.1) 
the classical stress energy tensor is 

Tab(X) 



y/tfx)Sg<*(x) 

= (1 - 2f) 4>- a 4>;b + - 4>-p4>'' P gab + 2£0 {4>- p P - 4>:ab9ab) 

+4? t [Rub - ^Rgab^j - \m 2 4? g ab (3.2) 

When we make the transition to quantum field theory, we promote the field 4>{ x ) to a field operator 4>{x). The 
fundamental problem of defining a quantum operator for the stress tensor is immediately visible: the field operator 
appears quadratically. Since 4>(x) is an operator-valued distribution, products at a single point are not well-defined. 
But if the product is point separated {4> 2 {x) — > 4>(x)4>(x')), they are finite and well-defined. 

Let us first seek a point-separated extension of these classical quantities and then consider the quantum field 
operators. Point separation is symmetrically extended to products of covariant derivatives of the field according to 

faa) 0;&) -» \ (g/Vp'Vb + g b p 'VaV P <) cj)[x)cj)(x'), (3.3) 
(<P;ab) -» \ (v Q V h + g* V'Vp'V,/) 4>(x)4>{x'). (3.4) 

The bi-vector of parallel displacement g a a (x,x') is included so that we may have objects that are rank 2 tensors at 
x and scalars at x' . 



To carry out point separation on (3.2), we first define the differential operator 



T a b = \ (1 - 2C) ( 5 a a 'V Q ,Vfc + g b b ' V a V^) + ^ - ^ gatg^'VcVd, 



6 The final expressions reaching such a length defy direct inspection check. However, since the algorithms are recursive, 
agreement at the lower orders extends to higher orders for the correctness of the complete result. For insurance, some of the 
larger expressions are derived by two independent methods and verified to be the same. 
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-£ (v a V b + g a a 'g b b ' V„/ Vy) + £.9a 6 (v c V c + V C ,V C ') 

^-Rab - ^.9ab^ - ^m 2 g ab (3.5) 

from which we obtain the classical stress tensor as 

T ab (x) = lim T ab <j){x)4>{x'). (3.6) 

x' —*x 

That the classical tensor field no longer appears as a product of scalar fields at a single point allows a smooth 
transition to the quantum tensor field. From the viewpoint of the stress tensor, the separation of points is an artificial 
construct so when promoting the classical field to a quantum one, neither point should be favored. The product of 
field configurations is taken to be the symmetrized operator product, denoted by curly brackets: 

m<t>{y) - \ {fa), kv)} = \ (4>{x)$(y) + Uy)${xj) (3.7) 
With this, the point separated stress energy tensor operator is defined as 

f ab (x,x') = ^T ab \j>{x)A{x')}- (3.8) 

While the classical stress tensor was defined at the coincidence limit x' — ► x, we cannot attach any physical meaning 
to the quantum stress tensor at one point until the issue of regularization is dealt with, which will happen in the next 
section. For now, we will maintain point separation so as to have a mathematically meaningful operator. 

The expectation value of the point-separated stress tensor can now be taken. This amounts to replacing the field 
operators by their expectation value, which is given by the Hadamard (or Schwinger) function 

GW(x,x') = ({4>(x),Hx')})- (3-9) 
and the point-separated stress tensor is defined as 

(f ab (x,x')) = l -T ab G^(x,x') (3.10) 

where, since T ab is a differential operator, it can be taken "outside" the expectation value. The expectation value of 
the point-separated quantum stress tensor for a free, massless (to = 0) conformally coupled (£ = 1/6) scalar field on 
a four dimension spacetime with scalar curvature R is 

(f ab (x, x')) = l - (gP' b G^-p'a + gP' a G« ;p , b ) - ^f' q G* 1 V g ab 

~ {9 P 'a9 q \G^, p , ql +G« ;a6 ) + 1 {(G^.y g ab ) 

+ I2 G(1) ( Rab -l R 9ab) (3-11) 



A. Finiteness of Noise Kernel 

We now turn our attention to the noise kernel. First introduced in the context of stochastic semiclassical gravity it 
is the symmetrized product of the (mean subtracted) stress tensor operator: 

8N ab , c , d ,(x,y) = {{f ab {x) - (f ab (x)),f c , d ,(y) - {f c , d ,(y))}) 

= ({f ab (x),f c , d ,(y)}) - 2(f ab (x))(f c , d ,(y)) (3.12) 

Since T ab {x) defined at one point can be ill-behaved as it is generally divergent, one can question the soundness of 
these quantities. But as will be shown later, the noise kernel is finite for y ^ x. All field operator products present 
in the first expectation value that could be divergent are canceled by similar products in the second term. We will 
replace each of the stress tensor operators in the above expression for the noise kernel by their point separated versions, 
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effectively separating the two points (x,y) into the four points (x, x',y, y'). This will allow us to express the noise 
kernel in terms of a pair of differential operators acting on a combination of four and two point functions. Wick's 
theorem will allow the four point functions to be re-expressed in terms of two point functions. From this we see that 
all possible divergences for y ^ x will cancel. When the coincident limit is taken divergences do occur. The above 
procedure will allow us to isolate the divergences and obtain a finite result. 

Taking the point-separated quantities as more basic, on e should replace each of the stress tensor operators in the 
above with the corresponding point separated version ( |3.8[ ), with T a b acting at x and x' and T c >d' acting at y and y' . 
In this framework the noise kernel is defined as 



8N ab , c < d <(x,y) = lim lim T ab T c , d < G(x, x' , y, y') 

x — >x y —*y 



(3.13) 



where the four point function is 

G(x,x',y,y') 



{{{4>(x),4>(x')},{kv),ky')}}) 
-2({kx),Hx')})({ky)J(y')}) 



(3.14) 



We assume the pairs (x, x') and (y, y') are each within their respective Riemann normal coordinate neighborhoods so 
as to avoid problems that possible geodesic caustics might be present. When we later turn our attention to computing 
the limit y — > x, after issues of regularization are addressed, we will want to assume all four points are within the 
same Riemann normal coordinate neighborhood. 

Wick's theorem, for the case of free fields which we are considering, gives the simple product four point function in 
terms of a sum of products of Wightman functions (we use the shorthand notation G xy = G + (x, y) = (<f>( x ) <M?/))) : 



{fa) 4>{v) H*') 4>(y')) = G xyl G yx > + g xx , g v 



G X y G x 'y' 



(3.15) 



Expanding out the anti-commutators in (3.14) and applying Wick's theorem, the four point function becomes 



G(x, x , y, y ) — G xy > G x 'y + G xy G x f y f 



Gy X ' Gy'x 



Gyx Gy'x' 



(3.16) 



We can now easily see that the noise kernel defined via this function is indeed well defined for the limit (a/, y') — + (x, y): 

G(x,x,y,y) = 2 {G 2 xy + G 2 yx ) , (3.17) 

^From th is we can see that the noise kernel is also well defined for y ^ x; any divergence present in the first expectation 
value of ( 3.14 ) have been cancelled by those present in the pair of Green functions in the second term. 



B. Explicit Form of the Noise Kernel 



We will let the points separated for a while so we can keep tract of which covariant derivative acts on which 
arguments of which Wightman function. As an example (the complete calc ulati on is quite long), consider the result 
of the first set of covariant derivative operators in the differential operator (3.5), from both T a b and T C 'd', acting on 
G(x,x',y,y'): 



i (1 - 20 2 (g a p "v p „ V h + g b p "v p „V a ] 

x (.g c < 9 "V^V rf , +gd> q '"v q i>>Vc>) G{x,x',y,y') 



(3.18) 

(Our notation is that V a acts at x, V c ' at y, Vh" at x' and V<j"' at y'). Expanding out the differential operator 
above, we can determine which derivatives act on which Wightman function: 

(1-20 2 



9c p '" 


g q " 


a {G xy ' | bp'" 


Gx'y;q"d' ~\~ G X y\bd' G x 'y'- 


q"p"' 








H~ Gyx' ; q 


' d' Gy' X -bp"' H~ Gyx-bd' Gy 


x'-q"p' 







"g q 


a {G X y' ; bp" 


G x' y.q" c' Gxy.bc' G x 'y' 


q"p"' 








H~ Gyx' ; q 


' C l Gy> X ]bp'" H" Gy X \bc' G y' 


x' \q" p" 


') 


+gc p ' 


' g q 


b {Gxy' ; ap" 


G x' y.q" d' H~ G X y\ad' G x 'y' 


■q"p"' 
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~h Gy X ' \q" d' Gy' X ;ap'" "i~ Gy X ;ad' G y' x' ;q" p'") 
~\~9d'^ 9^ b (G X y' \ap"' Gx'y-q" c' H~ G xy -ac' G x 'y' ;q" p'" 

~\~ Gyx' -q" c' Gy' x;ap"' ~t~ Gy X -ac' Gy' x ' \q" p"'Y\ (3.19) 

If we now let x 1 — > x and y' — > y the contribution to the noise kernel is (including the factor of i present in the 
definition of the noise kernel): 

g ^0 {Gxy.ad' G X y;bc' "h G X y,ac' G X y,bd') 

+ (1 — 2 £) {Gyx-ad' Gyx-bc' + Gyx-ac' Gyx-bd')^ (3.20) 

That this term can be written as the sum of a part involving G xy and one involving G yx is a general property of the 
entire noise kernel. It thus takes the form 

N aWd ,(x,y) = N aWd > [G+{x,y)] + N aWdl [G+(y,x)] . (3.21) 

We will present the form of the functional N abc i d i [G] shortly. First we note, for x and y time-like separated, the above 
split of the noise kernel allows us to express it in terms of the Feynman (time ordered) Green function Gf(x,v) and 
the Dyson (anti-time ordered) Green function Go{x,y): 

Nabc>d>(x,y) = N abdd , [G F (x,y)]+N abc , d , [G D {x,y)\ (3.22) 

7 The complete form of the functional N a bc'd> [G] is 

N a bc>d' [G] = N abdd , [G] + g ab N c>d , [G] + g c > d >N' ab [G] + g ab g dd >N [G] (3.24a) 



with 



&N abc id> [G] =(1 — 2£) (G [C ' b G- d 'a + G-c'a G- d ' b ) + 4£ 2 (G- C ' d ' G- a b + G G- abc 'd') 
— 2£ (1 — 2£) (G- b G- C 'ad' + G.a G :c i b d> + G ]d > G-. abc < + G- c > G- a bd') 
+2 £ (1 - 2 (G, a G. b R &d , + G, d G, d , R ab ) 

-4£ 2 (G, ab R c , d '+G, c/d ,R ab )G + 2eRc' d 'RabG 2 (3.24b) 
8N' ab [G] =2 (1-20 ((^-5) G, plb G, p ' a + ^{G., b G, p ,/ + G, a G, pl /) 



-h( ( 2 ^~l) G - p ' G -,*bp> +t; (c, p y G, ab +GG. abp y 

-(m 2 +£R') ((1-2 G- a G- b — 2 G ^ G-ab) 



+2£ ( (2£--J G.yG. p + 2G£G., p , p ) R ab 



1 ' £#) ZRabG 2 (3.24c) 



87V 



[G] =2(2£-~J G, plq G^ + Ae (G.y G ;q " + GG., p %y 



7 This can be connected with the zeta function approach to this problem |J as follows: Recall when the quantum stress tensor 
fluctuations determined in the Euclidean section is analytically continued back to Lorentzian signature (r — > it), the time 
ordered product results. On the other hand, if the continuation is r — > —it, the anti-time ordered product results. With this 
in mind, the noise kernel is seen to be related to the quantum stress tensor fluctuations derived via the effective action as 

= AT 2 abc , d , | t= _ ir>t , = _ iT , + ATt bc , d , \ m , =irl . (3.23) 
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(2£ - ((m 2 + £R) G., p , G- p + (m 2 + £R') G ;p G ;p ) 
! £ ((to 2 + C-R) G ;P ' P ' + (m 2 + G ; /) 



G 

i (to 2 + f J2) (m 2 + £R') G 2 (3.24d) 
For a massless, conformally coupled scalar field (£ = g an d m = 0), the noise kernel functional is 
72N a b C 'd' [G] = 4 (G ;c 'b G ;c i' a + G ;c ' a G-d'b) + G ;c /d< G ;a {, + GG- a bc'd' 

— 2 (G ; b G-c'ad' + G ;a G-c'bd' + G ; rf' G-.abd + G ;c ' G- a bd') 

+ 2 (G ;a G;fc i? C 'rf' + G; C ' G;<j' -R fl b ) 

— (G ;a 6 i? c 'd' + G ;c /(2/ i? a fc) G — -R c 'd' R a bG 2 (3.25a) 

288N' ab [G] — 8 ^—G-p'bG- p a + G-bG-p' a p + G- a G- P 'b p ^j 

4 (^G- P G- a bp' — G- p ' P G-ab + GG- a bp' P J 
— 2R' (2 G-a G-b — G G-ab) 

-2 (G, pl G p ' ~2GG, P A Rab-R'RabG 2 (3.25b) 



288iV[G] = 2G ;p - g G ; ;p ' <? + 4 (G^^G^' + GGp 
-4 (G ;p G ; ^ 9 ' + G/ G ;9 % 



P ,9 
9 



+i?G ;p , G ; p + i?'G ;p G ;p 

-2 (r G, p ' p ' + R' G ; /) G +\ R r ' g2 ( 3 - 25c ) 
For the minimal coupling (£ = 0) case: 

8N a bc'd' [G] = G ;C 'b G ; d'a + G ;C ' a G ;( j'b (3.26a) 

8N ab [G] = -Gyb G p 'a - m 2 G, a G, b (3.26b) 

87V [G] = i (G ;p / 9 G ; p ' 9 + 77i 2 (g ;p - G ; p ' + G ;p G ;p ^ + m 4 G 2 ) (3.26c) 



C. Trace of the Noise Kernel 

One of the most interesting and surprising results to come out of the investigations undertaken in the 1970's of the 
quantum stress tensor was the discovery of the trace anomaly ||(| . When the trace of the stress tensor T — g ab T a b is 
evaluated for a field configuration that satisties the field equation 

m = O^C-Ci?-- 2 )0 = O, (3.27) 

the trace is seen to vanish for the massless conformally coupled case. When this analysis is carried over to the 
renormalized expectation value of the quantum stress tensor, the trace no longer vanishes. Wald |2i| showed this was 
due to the failure of the renormalized Hadamard function G Ten (x, x') to be symmetric in x and x' , implying it does 



not necessarily satisfy the field equation (3.27) in the variable x' . (We discuss in the next section the definition of 
G ren (x,x') in the context of point separation.) 

With this in mind, we can now determine the noise associated with the trace. Taking the trace at both points x 



and y of the noise kernel functional (3.22) 
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N[G] =g ab g c ' d ' N abc , d , [G] 

= -3 G£ | (m 2 + ±£Rj G.y + (m 2 + ^R 1 ) G,/ 



~- G,/ + GG//J+ (m 2 + ItR) (m 2 + ^R'j G 2 

_3 (I ~ e ) ( +3 Q " e ) G;p> gpp ■ 3? ( G;p G;p,pp ' + Gp ' G/p ) 



i 2 + l^Rj G :p ,G/ +fm 2 + ^R>) G :i ,G l "\ (3.2s, 



For the massless conformal case, this reduces to 



N [G] = — {i?i?'G 2 - 6G (JZD ; + i?'D) G + 18 ((□(?) (D'G) + □'□(?)} (3.29) 
which holds for any function G(x,y). For G being the Green function, it satisfies the field equation ( 3.27 ): 



□ G=(m 2 +eR)G (3.30) 

We will only assume the Green function satisfies the field equation in its first variable. Using the fact D'R = 
(because the covariant derivatives act at a different point than at which R is supported), it follows that 

□ '□G= {m 2 +£,R)D'G. (3.31) 

With these results, the noise kernel trace becomes 
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x {g 2 (2m 2 + R'£) + (1-60 G. tP > G, p ' - 6 G£ G-y p ' } 
A {3 (2 to 2 + J? f ) G. p G ;p - 18 C G ;p G ;p ^ p ' 



in (i-c) G ;/; ,r;.'""[> c!.:;2) 



It vanishes for the massless conformal case. We have thus shown, based solely on the definition of the point separated 
noise kernel, there is no noise associated with the trace anomaly. Our result is completely general since we have 
assumed the Green function is only satisfying the field equations in its first variable. This condition holds not just 
for the classical field case, but also for the regularized quantum case, where we do not expect the Green function to 
satisfy the field equation in both variables. One can see this result from a simple observation: Since the trace anomaly 
is known to be locally determined and quantum state independent, whereas the noise present in the quantum field is 
non-local, it is hard to find a noise associated with it. This is in agreement with previous findings derived 
from the Feynman- Vernon influence functional formalism |l5[ . 

IV. REGULARIZATION OF THE NOISE KERNEL 

We pointed out earlier that field quantities defined at two separated points possess important information which 
could be the starting point for probes into possible extended structures of spacetime. Moving in the other (homeward) 
direction, it is of interest to see how fluctuations of energy momentum (loosely, noise) would enter in the ordinary 
(point-wise) quantum field theory in helping us to address a new set of issues such as a) whether the fluctuations to 
mean ratio can act as a criterion for the validity of semiclassical (note it is not that simple, see [Q), b) whether the 
fluctuations in the vacuum energy density which drives inflationary violates the positive energy condition, c) deriving 
structure formation from quantum fluctuations, or d) deriving general relativity as a low energy effective theory in 
the geometro- hydrodynamic limit ||. For these inquires we need to examine the consequences of taking or reaching 
the coincident limit and to construct regularization procedures to treat the ultraviolet divergences in order to obtain 
a finite result for the noise kernel in this limit. 

We can see from the point separated form of the stress tensor ( 3. 10] ) what we need to regularize is the Green function 



). Once the Green function has been regularized such that it is smooth and has a well defined x' — > x limit, 
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the stress tensor will be well defined. In Minkowski space, this issue is easily resolved by a "normal ordering" 
prescription, which hinges on the existence of a unique vacuum. Unfortunately, for a general curved spacetime, there 
is no unique vacuum, and hence, no unique mode expansion on which to build a normal ordering prescription. But 
we can still ask if there is a way to determine a contribution we can subtract to yield a unique quantum stress tensor. 
Here we follow the prescription of Wald [^6|, and Adler et. al. [^7j (with corrections ^8|) summarized in ||. We give 
a short synopsis below as it will be referred to in subsequent discussions and later papers. 
The idea builds on the fact that for G(x,x') LJ ~ (uj\cf>(x)(f>(x')\uj} , the function 

F(x, x 1 ) = G(x, x') Ul - G[x, x% 2 (4.1) 

is a smooth function of x and x' , where w\ and UJ2 denote two different states. This means the difference between the 
stress tensor for two states is well defined for the point separation scheme, i.e., 

Tab = I Hm T ab (F(x, x') + F(x', x)) (4.2) 

Z x' — >x 

is well defined. So a bi-distribution G L (x,x') might be useful for the vacuum subtraction. At first, it would seem 
unlikely we could find such a unique bi-distribution. Wald found that with the introduction of four axioms for the 
regularized stress tensor 

(f ab (x)) Icn = lim \% b (gV(x,x') - G L (x,x')) (4.3) 

x — >x Z V / 

G L (x,x') is uniquely determined, up to a local conserved curvature term. The Wald axioms are |2q,p|]: 



1. The difference between the stress tensor for two states should agree with ( |4.2| ); 

2. The stress tensor should be local with respect to the state of the field; 

3. For all states, the stress tensor is conserved: V a (T Q &) = 0; 

4. In Minkowski space, the result (0|T a b|0) = is recovered. 

We are still left with the problem of determining the form of such a bi-distribution. Hadamard |37| showed that 
the Green function for a large class of states takes the form (in four spacetime dimensions) 

G L (x, x') = -^( 2U ^ X "> + V (x, x') loga + W(x, x')) (4.4) 
87r z \ a J 

with U(x,x'), V(x,x') and W(x,x') being smooth functions. ^ We refer to Eqn (4.4) as the "Hadamard ansatz". 
Since the functions V(x,x') and W(x,x') are smooth functions, they can be expanded as 

00 

V(x, x') = v n (x, x')a n (4.5a) 

n=0 



W(x, x 1 ) = w n{x, x')cr n (4.5b) 

n=0 



with the Un's and w n 's themselves smooth functions. These functions and U(x,x') are determined by substituting 
G L (x, x') in the wave equation KG L (x, x 1 ) — and equating to zero the coefficients of the explicitly appearing powers 
of a n and <r™ log a. Doing so, we get the infinite set of equations 



8 When working in the Lorentz sector of a field theory, i.e., when the metric signature is ( — ,+,+,+), as opposed to the 
Euclidean sector with the signature (+,+,+,+), we must modify the above function to account for null geodesies, since 
a(x,x') — for null separated x and x' . This problem can be overcome by using a — > a + 2ie(t — t') + e 2 . Here, we will work 
only with geometries that possess Euclidean sectors and carry out our analysis with Riemannian geometries and only at the 
end continue back to the Lorentzian geometry. Nonetheless, this presents no difficulty. At any point in the analysis the above 
replacement for a can be performed. 
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with 



U{x,x') = A 1 / 2 ; (4.6a) 

2H v + KA 1/2 = 0; (4.6b) 

2nH n v n + Kv n -i = 0, n > 1; (4.6c) 

2H 2n v n + 2nH n w n + Kw n -i = 0, n > 1 (4.6d) 



//,. \ ( n - 1 - 1 ; (4.(>. ; ) 



^From Eqs (4.6d), the functions v n are completely determined. In fact, they are symmetric functions, and hence 
V(x, x') is a symmetric function of x and x 1 . On the other hand, the field equations only determine w n , n > 1, leaving 
Wo(x,x') completely arbitrary. This reflects the state dependence of the Hadamard form above. Moreover, even if 
wq(x 7 x') is chosen to be symmetric, this does not guarantee that W(x, x') will be. By using axiom (4) wq(x, x') = 0. 
With this choice, the Minkowski spacetime limit is 

(2ir) z a 

where now 2a = (t — t 1 ) 2 — (x — x') 2 and this corresponds to the correct vacuum contribution that needs to be 
subtracted. 

With this choice though, we are left with a G (x,x r ) which is not symmetric and hence does not satisfy the field 
equation at x' , for fixed x. Wald |^8| showed this in turn implies axiom (3) is not satisfied. He resolved this problem 
by adding to the regularized stress tensor a term which cancels that which breaks the conservation of the old stress 
tensor: 

PT) = (TT) + y^9* [«i] (4.8) 



where [v\\ = V\{x,x) is the coincident limit of the n = 1 solution of Eq (4.6d). This yields a stress tensor which 
satisfies all four axioms and produces the well known trace anomaly (T a a ) — [v\\ /8tt 2 . We can view this redefinition 
as taking place at the level of the stress tensor operator via 

Tab -> f ab + 37^35 9ab M (4.9) 

Since this amounts to a constant shift of the stress tensor operator, it will have no effect on the noise kernel or 
fluctuations, as they are the variance about the mean. This is further supported by the fact that there is no noise 
associated with the trace. Since this result was derived by only assuming that the Green function satisfies the fiel d 
equation in one of its variables, it is independent of the issue of the lack of symmetry in the Hadam ard a nsatz ( |4.4| ) . 



Using the above formalism we now derive the coincident limit expression for the noise kernel (3.24). To get a 
meaningful result, we must work with the regularization of the Wightman function, obtained by following the same 
procedure outlined above for the Hadamard function: 

G, cn (x,y) = G Ien , + (x,y) = G+{x,y) - G L (x,y) (4.10) 

In doing this, we assume the singular structure of the Wightman function is the same as that for the Hadamard 
function. In all applications, this is indeed the case. Moreover, when we compute the coincident limit of N a }, c id', 
we will be working in the Euclidean section where there is no issue of operator ordering. For now we only consider 
spacetimes with no time dependence present in the final coincident limit result, so there is also no issue of Wick 
rotation back to a Minkowski signature. If this was the case, then care must be taken as to whether we are considering 
[Nabc'd' [G ro „,+(x,y)]] or [N abc , d i [G Ien>+ (y, x)]]. 



We now have all the information we need to compute the coincident limit of the noise kernel (3.24). Since the point 
separated noise kernel N a t )C i a ''{x,y) involves covariant d eriv atives at the two points at which it has support, when we 
take the coincident limit we can use Synge's theorem ( [2.3| ) to move the derivatives acting at y to ones acting at x. 
Due to the long length of the noise kernel expression, we will only give an example by examining a single term. 



Consider a typical term from the noise kernel functional (3.24) 



Gyqyi\c' b G Ycn ;d' a ~f~ G rcn - C ' a G Tcn -d'b (4T1) 
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Recall the noise kernel itself is related to the noise kernel functional via 

Nabdd' = N abc , d > [G ren (x, y)] + N abc , d , [G ren (y, x)\ . (4.12) 

This is implemented on our typical term by adding to it the same term, but now with the roles of x and y reversed, 
so we have to consider 

Gren; c' b G rcn ; d' a "t - G r en • a > d G rcn ; b' c ~t~ G r en; c' a G rcn • d' b ~t~ Grcn;a' c G rcn ; b' d (4. 13) 

It is to this form that we can take the coincident limit: 

[G r en;c'6] [Cren;d'a] ~l~ [G r en;a'd] [G r en;&'c] ~\~ [G r en;c'a ] [G rcn -d'b] + [G ren; a' c ] [G ren ., b 'd} (4.14) 

We can now apply Synge's theorem: 

([G rcn;a];d [G ren ;ad]) ([G re n;h];c [G ren ;fcc]) 

[G 

rcn ; a d]) ([G rC n:c];6 — [G rcn ;6c]) 
+ ([G rcn ; a J ; c L^rer^a c]) ([G T en;b];d — [G rcn; bd\ ) 
H~ ([Crcn;c];a 

[G ron;QC ]) ([G ren ;d];b — [G re n;bd]) ■ (4-15) 

This is the desired form for once we have an end point expansion of G ron , it will be straightforward to compute the 
above expression. We will discuss the details of such an evaluation in the context of symbolic computations in a 
companion paper in this series. 

The final result for the coincident limit of the noise kernel is broken down into a rank four and rank two tensor and 
a scalar according to 

[Nabc>d>] = Kbcd + 9abK'd + 9cdKb + 9ab9cdN"' (4.16a) 

The two tensors and the scalar have the following form in terms of the concident limit and derivatives of the renor- 
malized Green function and its derivatives: 

8N a i, c d = (1 — 2£) {([G ren;a ]- C ~ [G ren -ac]) ([G r en;b];d ~ [G r en:bd]) 

c]) ([G r en;d];b ~ [G r en;bd]) 
-\- {[Gren;a];d ~ [G r en:ad]) ([Cren;6];c [G ren ;b c ]) 
~h ([GVen;d]:a — [Gren;ad]) ( [G ren ■ c ] : 5 — [G ren ■ b c ] ) } 

+4£ 2 {(-[G 

ren ; cd\ 

ren;abc\\d y^ren^a bd];c + [G ren;ab\\cd 

+ [G 

ren ; a bde]) 

~\~[G ren ] (— [G ren ;cda] ;b ~~ [G ren - C db\]a H~ [G ren - jC d\;ab H~ [G ren - : 

cdba})} 

~ 2£ (1 — 2 £) {([G ren ] ;a — [G ren;a ]) ([G ren:c d]-b ~ [Gren:cdb]) 

+ ([G ren ];b ~ [G r en:b\) ([G ren :cd\:a [G ren : cda\) 
H~ ([^Ven];c [Cren;c]) ([^ren;a6];d ~~ [G ren • abd\ ) 
H~ ([Cren];cZ ~~ [Gren;d\) ([^Ven;a&];c — [^Ven;a6c]) 
~f~[GVen;a] ( [^Ven;6c];<i [^ren;6d];c ~l~ [^ren;6];cd ~l~ [Cren;6dc]) 
H~[^ren;&] {~[G r en:ac\:d ~ [G*ren;a d]:c + [G ren ; a J ; cd + [Gren-.adc]) 
d];b - [G r en;bd\;a H~ [^ren;d];a& H~~ [^-?rero;6c2a]) 
H~ [^-* ren ; d] ( [^ren;oc] ; 6 [^Ven;6c];a [G ren:c \;ab H~ [^ren:6ca])} 
+2e(l-20{(([G re „] ;a [GVen;a]) ([^?ren];6 [Cren;6]) H~ [G r en;a ] [G re n;b]) Red 

+4C 2 {([G 

ren;a];6 H - [Cren;6];a [^?ren];a6 2 \Gren\a b }) Red 
H~ ([^ren;c];d H" [C'ren ; <i] ; c ~~ [G r e7i\;cd 2 [C ren;cc /]) R a b 

+ RabRcd[Gren}}[Gren] (4.16b) 

8N ab = 2 (2^ - (1 - 20 (([G ren . b ]. p - [G rm , pb ]) i\G ren , a r - [G ren .,*]) 

+ ([G ren;p\;b y^ren-.p b)) ([G re „ ;p ] ;a - [G ren . a p ])) 

+2^ (1 — 20 (([Gren];a — [G ren;a ]) ( [G ren :p P ] ■ b ~ [G ren -p P b] ) 

+ ([Gren];b [Gren:fc]) ([G re n;p^];a [Gren;p^a]) 

+[G ren;a } HG 
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+ [G ren; b\ ( [Gren:pa]' P ~ [G r en:a P ]:p + [G r en ; a] -p 1 '+ [Gren:a P p])) 

-4£ ( 2 £-0([G re7 ^] (-[G 

ren;paj;b [^"ren ;p6j ; a + [G ren;p];afe "I - [G r en;p6a]) 

+ ([G re „P- [G re „ ;J >]) ([G 

ren ; afr] ;p [G r en ; a bp])) 
—2 (rn 4- £ [G re n] ([G re n;a];& + [G ren; ti] ;a — [G ren ] ;a & — 2 [G re n;a6]) 

H~2 ^ I 2 £ —J ( [Gren] ;p [Gren] 2 [Gren] ;p [Gren ,P ] "4~ 2 [Gren ;p] [Gren ' P J ) ^?afc 



2, 

-4 £ [Gren] (2 [Gren;p]'^ [G r en];p^ 2 [G r en;p^]) ^afc 

-2 (m 2 +eR) ^afc[G re „] 2 



r 2 



87V'" =2^--j (([G ren;p ] ;g - [G re „ ;p ,]) ([G ren ;p p - [G re „^]) 

H~ ([Gren;g];p [G r en;pg]) ([Gren'^]'^ [Gren'^])) 
"^^^ ((^ [Gren;p] ~ [G r en];p^ [G r en;p p ]) [G ren;q J 

[Gren] (2 [Gren;p^g]'^ — [G r en ;p^] ; t/^ — [G r en ;p^t/^] ) ) 

— 8 £ ^2 £ — ^ (([G r en];p [G r en;p]) [Gren;q^] — [G r en];p [G r en;q^] ^ 

~l~ [G r en ;p 

} (2[G 

ren\q > \ < ^ H~ [G r en;q^] [G r en '^];q^ [G r en;t/^])) 

— 2 (m 2 + £R) ^2 £ — — ^ ([Gren];p [G ren ]' P — 2 [G re n];p [G re n' P ] + 2 [G re n;p] [Gren 1 

+4 (m 2 + i (2 [G ren;P P - [G re „] ; / - 2 [G re „ ; /]) [G re „] 
+ (m 2 +^) 2 [G re „] 2 

For the conformal coupling (£ = g and m = 0) case: 

72 A ^ abed 4 { ( [G ren ■ a ] j c [G ren ; aC ]) ( [G ren ; fc] ; [G ren ;bd]) 

+ ([G ren ; c] ; a [G r en ;ac ]) ([G ren ;d];6 ~ [G r en;bd]) 
+ ([G ren ; a ];d - [G ren : ad\ ) ( [G r en ; b\ : c [G r en : bc\ ) 
+ ([G ren;d];a [G r en;acz]) ([G r en:c]:fc [G r en;6c])} 
+ {(-[G ren; a]; 6 [G r en;6];a [G r en];afc H~~ [G r en;a 

*]) [G ren : caj 

-|-[Gren;afc] ( [Gren;c];c2 [Gren;d];c H~ [^ren];cd [Gren]cd\) 
H~[Gren] ( — [G r en; abc\ ; d [G r en;a 6d];c + [G ren;aoJ;ca ~r [ l - r ren;a bde]) 
H~[Gren] ( [^^reri. ; crfa] ; 6 [Gren;cdb\;a H~ [G re7l j C d];a6 ~l~ [G r en;cdfca ])} 
— 2 {([G r enJ;a [G r en;a]) ([Gren;c<i];& [Gren;ce£b]) 
H~ ([G r en];b [G r en;fc]) ( [G r en ; cd] ; a [G r en;cda]) 
H~ ([G r en];c [G r en;c]) ( [G r en ;a&] ; d [G r en ; a&d] ) 
H - ([Gren];rf [Gren;<i]) ([G r en;a6];c [G r en;afcc]) 

-|- [Gren; a] ( [G r en;6c];d [G r en;fed];c H - [Gye^i^jjcci ~l~ [Gren;6dc]) 
H~[Gren;6] {~[G ren:a c\:d [G r en;ad];c ~^ [^ren;o];cd H" [G r en;a dc]) 

~t~[G ren ; C ] ( [G ren ; a ^];5 [G ren ;5^]; a "t - [G ren ;rf]; a fo -\~ [G ren ;{,da]) 

+[G ren;c ij (— [G 

ren;oc];6 [Gy e 7i,|5c];o ~l~ [Gyen;c];o6 ~l~ [Gren;fcca])} 

+2{(([G ren] ; a [G r en ; a ]) ([G ren ] ; 6 — \G ren -b\) + [G 

ren;a] [G r en;b]) Red 
(,(,[G ren \-,c [Gren;c]) ([G r en];d [G r en;d]) H~ [G r en;c] [G r en;d]) Rab} 
{(\G ren ^a\;b H~ [Gren;6];a [G r en];ab — 2 [G r en;afc]) -^cd 
H~ ([G r en;c];d "I - [^ren;d];c [G r en];cd 2 [G r en;crf]) -^afe 
+ -R a 6 ^ccz[Gren]} [G ren\ 

UAN' ab = -4 {{[G ren , b \, p - [G ren;ph ]) ([G rt 

+ ([G ren ;pj ; o L^ren ;p 



P - fG 

aj L'- T ren: 


a P ]) 


;p l - [G 

J;a L^ren: 


;a P ])) 


];6 [Gren;p J 


°6]) 



T ren;p 



yp _ [G ren -b P ] :p + [G r 

en;b];p^ H - [Gren;6^p]) 
H~[Gren:6] ( — [G r en;pa] [G r en:a^];p H~ [G r en;a];p^ H~ [G r en:a^p])) 

2([G re „-f] (-[G 

ren;pa];fc [G r en;pb];a H~ [G r en;p];afc + [G re n;pfca]) 
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H~ ( [Cren] [G ren '^J ) ( [G ren ; ab\ ;p \G ren ; a bp })) 
— (R [Gren\ {\G re n;a\:b H~~ [^?ren;fe];a [Cren];afe 2 [G ren ;&£>])) 

ren];p [G ren \^ 2 [G ren ];p [Cren'^] ~t~ 2 [G ren ;p] [G ren '^]) Rab) 
2 [G ren ] (2 [Gren;p]'^ [Cren];p^ 2 [Gren;p^]) -^ab 

-(i?i? Qb )[G re „] 2 (4.17b) 

2887V'" = 2 (([G re n;p];g — [G r en:pq]) ([G re „' P ]' 9 — [G ren ' P9 ]) 
-\- (\Gren;q\;p [Gren;pg]) ([Gren'^]' P [Gj-en'^j)) 

-8((2[G 

ren;p] [^ren] ;p^ [Cren ;p ren;g J 

H-8 (([Gren];p [Cren;p]) [Gj- e 7i;<y ] [Cren];p [^ren;§ ]' 

~f~[^ren;p] (2 [^ren;^]'^ [Gren;q^] ^ [^ren'^];g^ [G ren jg^])) 
2 7? ([G ren ]-p [G ren ]'^ 2 [G ren ]jp [G ren '^] 4- 2 [G ren ;p] [G^e^'^]) 

+4i? (2 [G ren . p \ p - [G r 

en];p^ 2 [(?ren;p p ])[G renj 

+i? 2 [G re „] 2 (4.17c) 



For the minimal coupling (£ = 0) case: 



( [G ren j a 


;c [G> en - ac ]) {\G ren - j0 


:d ~ [G ren; bd]) 


( [G ren • d 


a [G ren ; ac ]) ([Gren;d] 


6 — [G r en;bd]) 


( [G?-en ;a] 


\d [Cren;at/]) ([Gj-en;&J 


;c [Cren;fec]) 


([G ren -d] 


a [G ren j at ;]) ([G ren ; C ] 


; b [G ren -5 C ] ) 



(4.18a) 

^^ab — ([^rera;&];p [G ren ;pb]) ([G^e^a]'^ [G ren ; a ^]) 
([Gren;p];fc [G ren ;pb]) ( [G ren '^];a [G ren ; ) 
Tfi {\Gren\;a [G reri ; a ]) ([G ren ];b [C ren :fo]) 

-m 2 [G ren;a ] [G ren;b ] (4.18b) 



8./V — ^ ( ( [G ren ;p] ; g [G ren ;pg] ) ( [G ren ' ^ [G ren ) 



1 

2 

ren ; q] ;p [G ren ;pq] ) ( [G ren '^j [^ren'^^])) 
-fm ([G re n];p [G ren ]' P — 2 [G ren ] ; p [G re „ ;p ]+2[G 

ren;p] [G ren '^]) 
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^ 4 [G re „] 2 (4.18c) 



V. FURTHER DEVELOPMENTS 

In this article we have derived a general expression for the noise kernel, or the vacuum expectation value of the 
stress energy bi-tensor for a quantum scalar field in a general curved space time using the point separation method. 
The general form is expressed as products of covariant derivatives of the quantum field's Green function. It is finite 
when the noise kernel is evaluated for distinct pairs of points (and non-null points for a massless field). We also have 
shown the trace of the noise kernel vanishes, confirming there is no noise associated with the trace anomaly. This 
holds regardless of issues of regularization of the noise kernel. 

The noise kernel as a two point function of the stress energy tensor diverges as the pair of points are brought together, 
representing the "standard" ultraviolet divergence present in the quantum field theory. By using the modified point 
separation regularization method we render the field's Green function finite in the coincident limit. This in turn 
permits the derivation of the formal expression for the regularized coincident limit of the noise kernel. 

When the Green function is available in closed analytic form one can carry out an end point expansion according to 



(2.5), displaying the ultraviolet divergence. Subtraction of the Hadamard ansatz (4.4), expressed as a series expansion, 
(4.5), will render this Green function finite in the coincident limit. With this, one can calculate the noise kernel for a 
variety of spacetimes, as will be developed in Papers II, III. Here we give an outline of the program in this series. 

A common analytic approximation is the Gaussian |38| ], which, in ultrastatic spacetimes, provides a closed form 
expression for the Green function. In Paper II we make use of this form to evaluate the noise kernel for a massless 
scalar field in hot flat space and the Einstein universe. For hot flat space, the Gaussian Green function is exact, 
while for the Einstein universe, the exact Green function is known. We also work explicitly with the ultrastatic 
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metric conformal to the Schwarzschild metric. Though the approximate Green function is known to be a fairly good 
approximation for the stress tensor, i.e., to second order, we find the approximation breaks down at fourth order (the 
noise kernel needs up to four covariant derivatives of the Green function). This manifest in the noise kernel trace 
computed under Gaussian approximation failing to vanish; indeed the magnitude of noise kernel becomes comparable 
to the components themselves. This form of the Green function fails to satisfy the field equations at fourth order. 

In Paper III we take advantage of the simple conformal transformation property of the scalar field's Green function 
to compute the noise kernel for a thermal fields in a flat FRW Universe as it is obtainable from that of hot flat space 
by a time dependent scale factor, and likewise use the Einstein Universe result to obtain the noise kernel for a closed 
FRW Universe. The most interesting case is when we conformally transform to the Schwarzschild metric. Though 
our results based on the Gaussian approximation break down close to the event horizon, we can obtain reasonable 
results for the fluctuations of the stress tensor of the Hawking flux in the far field region and check with analytic 
results pg|] . Along the way, we verify our procedure by explicitly re-deriving the Page result |3j| for the stress tensor. 
We note that in Page's original work, the direct use of the conformal transformation was circumvented by "guessing" 
the solution to a functional differential equation. Our result is the first we know where the methodology of point 
separation was carried all the way through to the final result. That we get the known results is a check on our method 
and its correct implementation. 

We presently outline how the work presented in this article are used to compute the coincident limit of the noise 
kernel. In future articles in this series we will present the details and results. We assume we have an analytic closed 
form expression for the Green function; we focus on the Gaussian approximation |3S}| . As mentioned above, this Green 
function is expanded in an end point series and the Hadamard ansatz is subtracted, generating the series expansion 
of the renormalized Green function. At this point we could directly substitute this expansion in our expression for 
the coincident limit of the noise kernel ( 4.16| ). The resulting expression becomes quite large and it is hard to glean 
any physical meaning from the resulting expressions. Instead we let this be the point in the problem when an explicit 
metric is chosen. Once this is done, in the context of the symbolic computation, it is straightforward to determine the 
component values of the Green function expansion tensors. For d = 4, there is at most 70 unique values. From this 
we can readily generate all the needed component values of the coincident limits of the covariant derivatives of the 
Green function, along with t he co variant derivatives of the coincident limits. It is these explicitly evaluated tensors 
that are then substituted in ( Lit ) to get the final result. 

We must stress that though all this is done on a computer, no numerical approximations are used; all work is done 
symbolically in terms of the explicit functional form of the metric and the parameters of the field. The final results 
are exact to the extent that the analytic form of the Green function is exact. 

As will be seen in Paper II and III the basic procedures for generating the needed series expansions are recursive 
on the expansion order. For the noise kernel we need results up to fourth order in the separation distance. The well 
established work for the stress tensor is to second order. This provides a check of our code by verifying we always 
get the known results for the stress tensor expectation value. Once we know the second order recursion is correct, we 
know the algorithm is functioning as desired and the (new) fourth order terms are correct. This becomes particularly 
important when we consider metrics conformally related, as we get intermediate results of up to 1100 terms in length. 

When we are interested in spacetimes for which we know the Green function in a conformally related space time, 
we have modified the above procedure. The simple conformal transformation property of Green functions allows us to 
get the Green function in the physical metric we are interested in from that in the conformally related metric (for the 
Gaussian approximation, this is the conformally optical metric). The main obstacle to overcome is the subtraction 
of the Hadamard ansatz. The divergent Green function is defined in terms of the optical metric while the Hadamard 
ansatz in terms of the physical metric. We need to re-express the transformed optical metric in terms of the physical 
metric. The defining equations for the geometric objects (e.g. Eqns 2.4 and 2.8) on the optical metric are transformed 
to the physical metric and recursively solved. Now the Green function series expansion can be written solely in terms 
of the physical metric. It is at this point that the symbolic computation environment reigns; the fourth order term 
in the expansion of the renormalized Green function has over 1100 terms. On the other hand, it is exactly in this 
context that we re-derive the Page |5f| result for the stress tensor vacuum expectation value on the Scharzschild black 
hole. As we stressed above, since all the series expansions used are recursively derived, so with the lower order results 
checked, we know the resulting expression is correct. Now having the general expansion of the renormalized Green 
function, we proceed as we outlined above and choose a metric and compute the coincident limit of the noise kernel. 

Though our main focus in this article has been to lay out the groundwork for computing the coincident limit of the 
noise kernel, our result ( [3.24a| ) is completely general. This expression for the noise kernel only assumes a scalar field 
and can be used with or without first considering issues of renormalization of the Green function. Also, our result 
for the coincident limit, holds regardless of the choice of Green function or metric. It only requires that the 

Green function had a meaningful coincident limit. It can become the starting point for numerical work on the noise 
kernel, an avenue worthy of some deliberation. 
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